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MINOS is a long baseline neutrino oscillation experiment designed to make precision mea- 
surements of the neutrino mixing parameters associated with the atmospheric neutrino mass 
splitting. Using a neutrino beam from the Main Injector (NuMI) facility at Fermilab, it 
compares the neutrino energy spectrum for neutrino interactions observed in two large de- 
tectors located at Fermilab and in the Soudan mine in northern Minnesota at a distance 
of 735km. We describe initial results for electron neutrino appearance in MINOS after two 
years of data-taking. We observe 35 events in the Far Detector with a background predic- 
tion of 27 ± 5(stat.) ± 2(syst.) based on the measurement in the Near Detector. The 1.5cr 
excess of events can be interpreted in terms of ffi—^fe oscillations. At 90% CL we obtain 
an upper limit range of sin'^(26'i3) < 0.28 — 0.34 for the normal neutrino mass hierarchy and 
sin^(20i3) < 0.36 — 0.42 for the inverted hierarchy depending upon Sep- 

1 Introduction 

Over the last decade, several experiments have provided compelling evidence for neutrino oscilla- 
tions by detecting neutrinos originated in the atmosphereEEl or in the sunllH These observations 
have been confirmed more recently by a subsequent generation of experiments based on artificial 
sources of neutrinos, such as acceleratorsEEl and reactors^. The results support the description 
of three flavors of neutrinos changing identity as they travel and being related to three mass 
eigenstates as described by the Pontecorvo-Maki-Nakagawa-Sakata mixing matrix!^. 

The MINOS long baseline neutrino experiment at Fermilab has provided compelling evidence 
for disappearance, establishing that neutrinos oscillate as they travel^. MINOS has a baseline 
of 735 km and neutrinos peaking at 2-3 GeV, thus it is dominated by the atmospheric mass 
splitting of MINOS, in fact, provides the most precise measurement for this mass scale: 

|Am32| = 2.43ib0.13 x 10"'^ eV^. In this regime, the dominant oscillation channel is expected to 
be Vfj, — > i^T) but it is possible that — > i/g oscillations could occur. Observation of this oscillation 
channel would imply a non-zero value of the mixing parameter ^13, one of the missing pieces of 



the neutrino puzzle. A non-zero value of this parameter opens the door to the observation of 
CP violation of the neutrino sector. As of this writing, disappearance of the reactor Ve over a 
short baseline of few km has not been observed'^. This result implies that ^13 must be small, 
as measured by the Chooz experiment, sin^(26'i3) < 0.15 at 90% confidence level (CL) for the 
MINOS measured value of Amgg. Two other experiments have given limits, albeit with lower 
sensitivitylMm 

2 The MINOS Experiment 

In MINOS, neutrino interactions from the Fermilab Main Injector (NuMI) facility'^ are recorded 
at two detectors: Near Detector (ND), 1 km from the NuMI target, and a Far Detector (FD) 
at 735 km from the same target. The detectors are tracking calorimeters composed of 2.54 cm 
thick steel absorber planes and 1.0 cm thick active scintillator planes forming detector layers 
that correspond to ^ 1.4 radiation lengths. The scintillator planes are composed of 4.1 cm wide 
strips which corresponds to ~ 1.1 Moliere radii. The data recorded in the ND establishes the 
properties of the beam before oscillations have occurred^ Evidence for oscillations is observed 
as distortions of the beam spectrum or composition measured at the FD with respect to the ND. 
Thanks to the similarity between the two detectors, systematic errors arising from uncertainties 
in the neutrino interactions or in the neutrino flux largely cancel. 

The neutrino beam is produced by impinging 120 GeV protons from the Fermilab Main 
Injector upon a graphite target, producing pions and kaons, which are then focused by two 
magnetic horns The horn current and position of the target relative to the horns can be 
configured to produce different neutrino energy spectra. Most of the physics data has been 
acquired in the low energy configuration with a peak at 2-3 GeV, but other configurations have 
been used to study the backgrounds and systematics. The beam is 98.7% and I?^ and 1.3% 
intrinsic i/g+t'e- The latter originate from decays of muons produced in pion decays and from 
kaon decays. The muon component dominates below 8 GeV and it is well constrained by the 
z/^-CC flux from different beam configurations 0IH Errors on the fg content of the beam are 
3.5% at the ND and 5.7% at the FD. 

The FD data for the analysis presented here was recorded between May 2005 and July 2007, 
corresponding to an exposure of 3.14 x lO^'^ protons on target (POT). The ND data are sampled 
at a 1/6 rate uniformly over the 2 year run period. The data are partitioned in two subsets, with 
the second data set acquired with a ~ 1 cm longitudinal offset relative to the target position of 
the first data set. This offset results in a ~ 30 MeV shift in the neutrino spectrum that has been 
incorporated in the Monte Carlo (MC) simulation appropriately weighted by exposure. The 
MINOS MC simulation of the beam line and detectors is based on GEANTs'^. The simulation of 
the hadron production in the target is done using FLUKaIMH The NEUGEN simulation is used to 
model the neutrino interactions and it has been tuned to available bubble chamber data from 
relatively higher invariant masses than the region of interest for MINOS. The fragmentation 
or hadronic shower model uses the KNO description for the low invariant mass and transitions 
gradually at higher masses to the JETSET model13. 

3 Finding Electron Neutrinos in MINOS 

The flavor of the incoming neutrino can only be identified in charged current (CC) neutrino 
interactions, when the charged lepton partner of the incoming neutrino is produced in the final 
state. The MINOS detectors were designed to preferably detect long muon tracks resulting 
from muon- neutrino charged current interactions -|- N — > ^ + X (f^-CC). However, one can 
identify electron-neutrino charged current interactions i^e + N — > e + X (fe-CC) by searching 
for electrons which deposit their energy in a narrow and relatively short span in the MINOS 



calorimeters with a longitudinal distribution as described by the gamma function'^. Additional 
calorimeter activity can be produced by the hadronic showers resulting from the breakup of 
the recoil nucleus, X. There are other neutrino scattering processes that can produce similar 
topologies in the MINOS detectors, such as neutral current (NC) interactions v + 1^ ^ v + X , 
as well as i/^j-CC interactions where little energy is transferred to the out-going muon. Both 
of these background components are dominated by hadronic showers, typically possessing an 
electromagnetic element arising from vr'^ decays. Other less significant backgrounds, arise from 
beam v^-CC interactions, u-r-CC interactions from oscillations and cosmogenic sources. 

The data are first filtered to ensure data and beam quality. Fiducial volume cuts are imposed 
to assure well reconstructed events. In order to select a v^-CC like event sample, we restrict 
the energy range between 1 and 8 GeV, which contains the — > Ve oscillation probability 
maximum. The lower cut removes mainly NC events, while the higher cut removes beam v^-CC 
events. Background from cosmogenic sources are suppressed to less than 0.5 event (at 90% 
CL) by selecting events in time with the accelerator beam pulse and in the direction of the 
beam. Events are required to have a reconstructed shower and at least 5 contiguous planes with 
energy deposition greater than 0.5 MIP. Events with tracks longer than 25 planes or 12 track like 
planes are rejected. MC simulations indicate that these cuts improve signal to background ratio 
from 1:55 to 1:12 assuming a sin'^ (2^13) value at the CHOOZ limit. Two selection algorithms 
have been developed to enhance this ratio for higher background rejection, they are the ANN 
and LEM methods. The two selection algorithms rely on very different techniques, provide 
different signal to background ratios, and are sensitive to different systematic uncertainties. In 
the analysis reported here, the ANN selected sample is used to derive the final results but the 
LEM selection is presented as a cross check. 

3.1 The ANN selection method 

The final fg-selected sample is obtained using a method based on an Artificial Neural Network 
(ANN) with 11 variables characterizing the longitudinal and transverse energy deposition pro- 
files'^. Some of these quantities are the energy fraction in windows of 2, 4 or 6 planes, the 
fraction of energy in a 3 strip wide road, the RMS of the transverse energy deposition, etc. The 
neural network is trained to separate the signal Vg-CC events from NC or v^-CC background. 
The ANN generates a single numerical output with higher values indicating likely z^e-CC events. 
The acceptance cut value is chosen to be at 0.7 by maximizing the ratio of the accepted signal to 
the expected statistical and systematic uncertainty of the background. Assuming a signal at the 
CHOOZ limit, this method gives a 1:4 signal to background ratio, with a 41% signal efficiency 
and a rejection efficiency at 92.3% and 99.4% for NC and v^-CC respectively. 

3.2 The LEM selection method 

A second selection method, the Library Matching Method (LEM), is used as a cross check. In 
this novel technique, each candidate is compared to a large library of simulated Vg-CC and NC 
event^^. The 50 best matches are identified based on the relative probability that the hit pattern 
in a library and a candidate event come from the same neutrino interaction. Three variables are 
constructed: the fraction of these matches that are fg-CC events, the mean hadronic y of the best 
matches, and the mean fractional charge matched within those best matches. These quantities 
are then used as a function of energy to calculate a likelihood function that discriminates between 
signal and background. The acceptance cut value is chosen at 0.65. The LEM method gives 
better background rejection than the ANN, with a signal to background ratio of 1:3, but also 
displays increased sensitivity to certain systematic uncertainties. This method results in a 46% 
signal efficiency and a rejection efficiency at 92.9% and 99.3% for NC and t'^-CC respectively. 
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Figure 1: Reconstructed energy spectrum of z^e-selected ND data with total MC (left). The black points show data 
with statistical errors, the red histogram shows the MC with systematic errors as a shaded region. Reconstructed 
energy spectrum of z/g-selected ND data with predicted background components (right). The solid histogram 
corresponds to the total of the 3 background components present in the ND data and it has been tuned to agree 
with the data points. The selected background from NC is shown in blue and from v^-CC in red, both as obtained 
by the horn-off method. The shaded histogram shows the beam v^-CC component from the MC simulation. The 
errors in the data are statistical and are obscured by the markers; the errors on the components are systematic. 



4 Estimating the Background for Electron-Neutrino Appearance 

Data in the ND are recorded before oscihations have occurred, thus events selected there by 
the ANN are background events. The data recorded in the ND are directly used to predict 
the number of background events expected in the FD. The background measured in the ND 
is comprised of three different components: NC events, v^-CC events and beam i^g events. 
Oscillations and beam line geometry considerations require that each component be treated 
separately in the prediction of the FD backgrounds. Figure [l] (left) shows the energy spectrum 
of events selected with the ANN as z^g-CC like for data and MC in the ND. The data agrees 
with the simulation within the uncertainty which is dominated by the hadronic shower model 
systematic errors. 



4-1 Background estimation using horn-off data 

The background components are determined using a NC-enriched data sample recorded with 
the focusing horns turned off (horn-off). In this configuration the pions are not focused, the low 
energy peak of the neutrino energy distribution disappears, leaving an event sample dominated 
by higher energy NC events. These data used in conjunction with the standard beam configu- 
ration (horn-on) and the simulated horn-off to horn-on ratios allow us to extract the individual 
NC and u^-CC components as a function of energy. The beam fe-CC are used as an input for 
this method and are well constrained by the ND f^-CC data. The ratios are shown to be well 
modeled before the application of the z^g-CC selection and uncertainties such as the hadronic 
shower modeling cancel to first order in these ratios. Figure [l] (right) shows the NC, v^,-CC, 
and intrinsic beam z^g-CC components of the background in the ND derived used this technique. 
Systematic errors on the components arise from uncertainties in the beam flux, cross section 
and selection efficiency and are derived from the data for the main background components, NC 
and f^-CC. 



Table 1: Breakdown of systematic uncertainties contributed by different sources to the total number of background 

events in the Far Detector. 



Uncertainty source 


Uncertainty on 
background events 


Far/Near ratio (systematic) 


6.4% 


Horn-off (systematic) 


2.7% 


Horn-off (statistical) 


2.3% 


Total Systematic Uncertainty 


7.3% 


Expected Statistical Uncertainty 


19% 



Background estimation using muon-removed data 

A second technique allows us to study the ND measured background using an independent 
sample of showers derived from f^-CC events selected with long tracks'^. These events are 
processed to remove the hits associated with the muon tracklHHH Then the shower remnant is 
reconstructed to provide an independent sample of hadronic shower events. The procedure is 
done on data and MC and the Uf, selection is applied to both. The ratio of these muon-removed 
data to muon-removed MC is used to correct the NC component of the background. The residual 
difference between data and MC is absorbed into the f^-CC component by making the total MC 
agree with the data. The predictions for the background components from this method agree 
well with those obtained from the horn-off method. 

5 Far Detector Signal and Background Prediction 

Once the ND energy spectrum is decomposed into its background components, each of these 
spectra is multiplied by the Far to Near ratio from the MC simulation for that component 
providing a prediction of the FD spectrum. The MC simulation takes into account differences 
in the spectrum of events at the ND and FD due to beam line geometry as well as differences 
in detector calibration and topological response. Oscillations are included when predicting the 
v^-CC component. The smaller Vr-CC and beam fg-CC components are derived from the 
f^-CC selected FD event energy predicted spectrum. For the ANN selection, we expect 26.6 
background events, of which 18.2 are NC, 5.1 are f^-CC, 2.2 are beam z/g-CC and 1.1 are z^r-CC; 
using Am|2=2.43 x 10~^eV^, sin^ (2^23)=l-0, and sin^ (2^i3)=0 as oscillation parameters. 

It is possible to estimate the efficiency for selecting fg-CC events by using the sample of muon- 
removed events and embedding a simulated electron of the same momentum as the muon. Test 
beam measurements indicate that the ANN selection efficiency on single electrons agrees with the 
simulation to within 2.6% in the MINOS detectord^IHH Comparisons between the muon-removed 
(with added electron) data and MC simulation samples indicate that the selection efficiency of 
the signal events is well modeled. The fe-selection algorithms focus on the electromagnetic core 
of the shower and seem not to be affected by the hadronic shower. The difference in selection 
efficiency between data and MC is -0.3% and is used as a correction to the signal selection 
efficiency. 

The systematic errors were evaluated by generating modified MC samples and quantifying 
the change in the number of predicted background events in the FD^'^. Table [l] shows that the 
dominant uncertainty arises from Far /Near differences such as the relative energy scale calibra- 
tion, the details of the modeling of the photomultiplier tubes and the relative normalization. 
Other uncertainties arising from neutrino interaction physics, shower hadronization, intranu- 
clear re-scattering, and absolute energy scale errors largely cancel out in the extrapolation. The 
Far/Near systematic errors are added in quadrature along with the systematic and statistical 
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Figure 2: Distribution of the ANN selection variable for data and MC events in the FD (left). Black points show 
data with statistical error bars. The red histogram shows the background expectation. Below, the purple shaded 
histogram shows the z/e-CC prediction by the oscillation hypothesis compared to the data excess (black points). 
The reconstructed energy distribution of the v^-CC selected events in the FD (right). Black points show the 
data with the statistical errors bars. The stacked histogram shows the total predicted background, subdivided 
in components: NC (blue), t'^-CC (red), Ur-CC (green) and beam Vg-CC (magenta). Below, the purple shaded 
histogram shows Ve-CC prediction by the oscillation hypothesis compared to the data excess (black points). 

error arising from the background decomposition in the ND. As shown, the error is dominated 
by the expected statistical uncertainty on the number of background events for this data set. 

6 Far Detector Data 

We estabhshed the prediction of the background components and the signal as well as the 
systematic uncertainties before examining the data in the FD. In addition, we studied a signal 
free sample to validate the Vf, selection methods and background estimation techniques we have 
developed. All FD data was studied with both V(.-CC selection methods. 

6.1 Results using the ANN method 

The muon removal procedure was applied to the FD data and MC simulation. The number of 
events passing ANN selection was compared to the prediction based on the same sample from 
the ND. In this muon-removed sample, we observe 39 events above the selection cut, with an 
expectation of 29 it 5(stat.)ib2(syst.). The excess of events is concentrated at very high values 
of the selection variable. Examination of the events and distributions give no evidence that 
the excess is not a fluctuation. The sample will be explored further with the larger data set 
currently being acquired. 

The standard FD data sample was first studied for events passing all but the fg-CC selection 
criteria. In the FD data, 146 events were observed below ANN=0.55, compared to a pure 
background expectation of 132 it 12(stat.)ib8(syst.). In the signal region above ANN=0.7, we 
observe 35 events, with a background expectation of 27ib5(stat.)ib2(syst.). Figure [2] (left) shows 
the ANN selection variable distribution. The observed energy spectrum for events within the 
signal region is shown in Figure |2] (right). 




Figure 3: The 90% CL in the sin^ (2^13) and 5cp plane using the ANN selection method. Black lines show the 
best fit to the data for both the normal hierarchy (solid) and inverted hierarchy (dotted). Blue (red) lines show 
the 90% CL boundaries for the normal (inverted) hierarchies. 



6.2 Results using the LEM method 

The second selection method, LEM, was used as a cross check and it is not described in detail 
here. For the muon removed-sample using this selection method, we observe 25 events, with 
an expectation of 17 it 4(stat.)ib2(syst.). For the standard FD data in the region below a cut 
of LEM=0.55, we observe 176 events, with an expectation of 157 it 13(stat.)ib3(syst.). As in 
the ANN method, we observe a small excess (< 2a) for both the muon-removed sample and 
the region below the selection cut. In the signal region, we observed 28 events compared to 
a background expectation of 22 it 5(stat.)ib3(syst.). These results are also consistent with the 
ANN selection. 



7 Results for Ve in MINOS 

Figure [3] shows the 90% confidence intervals in the sin^ (2^13) and the CP phase, 5cp^ plane 
obtained using the observed total number of events selected by the ANN method. The oscil- 
lation probability is computed using a full 3-flavor neutrino mixing framework that includes 
matter effects'^ and introduces a dependency on the neutrino mass hierarchy, the sign of /S.m'^2- 
The contours are calculated using the MINOS best fit values of |Am|2| =2.43 x lO^^eV^ and 
sin^ (2^23) =1-0. Statistical (poisson) and systematic effects (gaussian) are incorporated via the 
Feldman-Cousins approach which we use to determine the desired confidence intervals'^Sl. 

In conclusion we have described the first results of a search for appearance in the MINOS 
experiment. The 1.5cr excess of events can be interpreted in terms of — > Vf, oscillations. At 
90% CL we obtain an upper limit range of sin'^(20i3) < 0.28 — 0.34 for the normal neutrino mass 
hierarchy and sin^(20i3) < 0.36 — 0.42 for the inverted hierarchy depending upon 5cp- 

MINOS is actively taking data and has already doubled the statistics used in this analysis 
(Spring 2009). In the near future expect an improved analysis taking advantage of the larger 
statistics (> 7 x lO^'^ POT ) and improvements in our understanding of the systematic errors. 
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